Spectral, optical and transport properties of the adiabatic anisotropic Holstein model: 
Application to slightly doped organic semiconductors 
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Spectral, optical and transport properties of an anisotropic three-dimensional Holstein model 
are studied within the adiabatic approximation. The parameter regime is appropriate for organic 
semiconductors used in single crystal based field effect transistors. Different approaches have been 
used to solve the model: self-consistent Born approximation valid for weak electron-phonon coupling, 
coherent potential approximation exact for infinite dimensions, and numerical diagonalization for 
finite lattices. With increasing temperature, the width of the spectral functions gets larger and larger 
making the approximation of quasi-particle less accurate. On the contrary, their peak positions 
are never strongly renormalized in comparison with the bare ones. As expected, the density of 
states is characterized by an exponential tail corresponding to localized states at low temperature. 
For weak electron-lattice coupling, the optical conductivity follows a Drude behavior, while, for 
intermediate electron-lattice coupling, a temperature dependent peak is present at low frequency. 
For high temperatures and low particle densities, the mobility always exhibits a power-law behavior 
as function of temperature. With decreasing the particle density, at low temperature, the mobility 
shows a transition from metallic to insulating behavior. Results are discussed in connection with 
available experimental data. 
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I. INTRODUCTION 

Organic field-effect transistors (OFETs) play a key role 
in the field of plastic electronics. Recently, single crys- 
tal OFETs have been developed showing charge mobil- 
ity at least one order of magnitude larger than that of 
thin films. ^ The most promising are those based on poly- 
acenes, such as pentacene and rubrene, that exhibit a 
strong anysotropy and the largest mobility measured in 
organic semiconductors.^ 

In spite of many applications based on such devices, 
the intrinsic transport mechanism acting in high mobil- 
ity OFETs is not fully understood. At low tempera- 
tures, an insulating behavior takes place with a sharp 
fall of resistivity. At temperatures higher than lOOif, 
the mobility of these systems exhibits a power-law 
behavior (^p ~ T~^)^— which cannot be simply as- 
cribed to band transporti^ Moreover, in some systems, 
starting from room temperature, there is a crossover 
from band-like to activated hopping behavior.®'^ Very 
recently, in order to better investigate the properties of 
the quasi-particles involved in the transport mechanism, 
many experiments based on spectroscopic probes such 
as photoemission^"— and optical conductivitj*i2rii have 
been performed in pentacene and rubrene. Photoemis- 
sion experiments suggest that the quasi-particle energy 
dispersion does not exhibit a strong mass renormalization 
even if the width of the peaks increases significantly with 
temperature. Optical conductivity measurements cannot 
be easily explained in terms of a Drude-like response, 
since they show a low frequency peak whose origin is not 
clear. Both kinds of experiments indicate moderate val- 
ues of electron-phonon (cl-ph) coupling pointing out that 
polaron formation should not play a prominent role for 
the interpretation of data. 

From the theoretical point of view, the calculated scat- 
tering lengths of the charge carriers are too short to 



be compatible with the band-like transport mechanism 
found at intermediate temperaturesi^ Therefore, in or- 
der to analyze the tunneling-hopping crossover in some 
transport data, a mechanism based on polaron band nar- 
rowing with increasing temperature has been proposedJ^ 
For the theoretical analysis, the simplest model for molec- 
ular crystals, the Holstein xnodel^t^ has been quite ex- 
tensively studied within several el-ph coupling regimes. 
Recently, in order to take into account the anisotropy 
of the electronic band and the most important couplings 
coming from the interaction of charge carriers with low- 
frequency phonons, a three-dimensional anisotropic Hol- 
stein model for a single electron has been analyzedJ^ 
However, a study of spectral and optical properties in 
the regime appropriate to organic semiconductors has not 
been performed up to now. 

For high mobility systems, a very interesting and sim- 
ple one-dimensional model has been recently introduced 
by Troisi and Orlandifi^ where the charge carrier in- 
teracts with low-frequency modes leading to a modula- 
tion of the charge carrier hopping. The model is pro- 
posed for an approximate description of rubrene and it 
is somehow very close to the Su-Schrieffer-Heeger (SSH) 
model introduced in a different context.'^*' The model has 
been studied within the adiabatic limit (phonons obey a 
classical dynamics) by using an approximated dynam- 
ical approach. Some results have been confirmed and 
extended by Ciuchi and Fratini,— who, instead, used a 
static approach, which is equivalent to the classical prob- 
lem of a quantum particle in the presence of disorder due 
to electron-lattice coupling. The main results of these 
studies is that a power law for the mobility temperature 
dependence can be recovered within the proposed model 
and that the charge carriers involved in the transport un- 
dergo a " dynamical localization" with increasing temper- 
ature. Very recently, some of us have made a systematic 
study of the transport properties of the SSH model in the 
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adiabatic limit by using the Kubo formula for the conduc- 
tivity showing that the inclusion of the vertex renormal- 
izations is relevant for the interpretation of experimental 
resultsj^ Moreover, in this paper, spectral and optical 
properties of finite density systems have been discussed 
in connection with experimental data. 

It is clear that one-dimensional adiabatic models suffer 
of severe limitations, of which the main is that electronic 
states are always localizedi^ Moreover, ingredients such 
as band anisotropy, small but finite carrier density are 
necessary for a correct description of the systems. There- 
fore, in this pap er, we study the Holstein model pro- 
posed in Ref.jlq within the adiabatic approach focus- 
ing on the weak to intermediate el-ph coupling regime 
which is relevant for high mobility organic semiconduc- 
tors. Different methods have been used to solve the 
model: self-consistent Born approximation (SCBA) accu- 
rate for weak el-ph coupling, coherent potential approx- 
imation (CPA) exact for infinite dimensions, and exact 
diagonalization (ED) for finite lattices. 

The calculated quantities are the spectral, optical and 
transport properties of the model. With increasing tem- 
perature, the width of the spectral functions gets larger 
and larger making the quasi-particles less defined. On 
the contrary, the spectral functions show peaks which are 
weakly renormalized in comparison with those of the bare 
bands. The marked width of the spectral functions gives 
rise to densities of states with a low energy exponential 
tail increasing with temperature. At low temperatures, 
this tail corresponds to localized states and gives rough 
indications for the energy position of the mobility edge. 
With increasing temperature, in the regime of low carrier 
doping investigated in this paper, the chemical potential 
always enters the energy region of the tail. These fea- 
tures allow to reconcile the results provided by the An- 
gle Resolved Photoemission Spectroscopy (ARPES) data 
(band-like description) with the trend that charge carri- 
ers appear more localized at high temperature. 

The crossover between low to high temperature be- 
havior considered in the spectral properties has a rele- 
vant effect also on the optical and transport properties. 
For weak el-ph coupling, the optical conductivity follows 
a Drude behavior and the mobility fully described by 
SCBA can be ascribed to band transport. In the in- 
termediate el-ph coupling regime, vertex corrections be- 
come important and an exact diagonalization is needed. 
In this regime, at low temperature, a temperature de- 
pendent low frequency peak appears in the conductiv- 
ity, while, at higher temperatures, the spectral features 
broaden over quite a large frequency range. The presence 
of this peak confirms the findings of a previous paper 
and compares well with recent experimental results fi^""— 
Finally, the mobility exhibits a power-law behavior for 
any high temperature. For very small particle density, 
at low temperature, the mobility shows a transition from 
metallic to insulating behavior with decreasing the par- 
ticle density. Results are discussed in connection with 
available experimental data. 



The paper is organized as follows. In section II the 
model and the methods of solution are presented, in sec- 
tion III spectral properties are analyzed, in section IV 
optical and transport properties are discussed. 

II. THE MODEL 

The spectral, optical and transport properties will be 
studied within the Holstein mode l^^'^^ in the adiabatic 
limit including the anisotropy important for the descrip- 
tion of organic semiconductors.^* It can be summarized 
in the following model hamiltonian: 

Ri Ri 

where and are the classical oscillator displace- 
ment and momentum, respectively, relative to the site 
Ri of a cubic lattice with parameter a, m the oscillator 
mass, k the elastic constant, and the electronic part Hei 
is 

Ri,S Ri 

In eq. ([2]), t^g^ is the bare electron hopping toward the 
nearest neighbors 6, ct, and cs are the charge carrier 

Ri ^ 

creation and annihilation operators, respectively, and a 
is the coupling constant that controls the link between 
the electron density and lattice displacement. We use 
units such that lattice parameter a = 1, Planck constant 
h = 1, Boltzmann constant fc^ = 1, and electron charge 
e = 1. 

The treatment is semiclassical since the electron dy- 
namics is fully quantum, while the ion dynamics is as- 
sumed classical. The latter approximation is well justi- 
fied from the typical values of phonon frequencies = 
•y/fc/m and hopping constants of high mobility organic 
semiconductors such as pentacene or rubrene. Following 
Ref.[ll], Wo ~ lOmel/, ~ lOOmeF, ~ bQmeV, ty ~ 
20meV leading to an adiabatic ratio 7 = wq/^z = 0.1. 
Therefore, all the results reported in this paper are valid 
for temperatures T > cjq — 12QK. 

In the adiabatic limit, the following quantity A 



correctly describes the strength of the el-ph coupling. Ac- 
tually, in the adiabatic limit, the calculation is equivalent 
to the classical problem of quantum particles in the pres- 
ence of diagonal disorder = axj^ due to electron- 
lattice coupling. The disorder is distributed according 
to the probability function of the lattice displacements 

P (^{xj^,}^ , that has to be self-consistently calculated as 

a function of electron-phonon coupling A, temperature 



3 



T and particle density c = Np/L, with Np number of 
particles and L number of sites of the cubic lattice. 

In this paper, we will focus on the weak to interme- 
diate electron-lattice regime that seems to be appropri- 
ate to systems like high mobility organic semiconduc- 
tors. Simulations will concentrate on temperatures lower 
than the typical electron energy of the order of (this 
temperature range is extensively investigated in many 
OFETs), and they will not explore the high tempera- 
ture regime.— 1^ Moreover, in most OFETs, the induced 
doping is not very high (typically much smaller than one 
charge carrier for ten molecules), therefore, in the fol- 
lowing, we will investigate the regime of low doping (up 
to c = 0.01). For this regime of parameters, the prob- 
ability function of the lattice displacements P ^{x^ }^ 

shows very tiny deviations from the distribution of free 
oscillators (it has been checked also in a recent paper^^ 
for a related model), so that one can correctly assume 



2m7r 



L/2 



exp 



-4e 



R^ 



(4) 



with /3 : 



1/T. The mean value of an observable 
over the distribution P {^{xj^,}^ is indicated 



by the standard symbol 



(5) 

In the next three subsections, the methods of solution 
for the considered model are discussed. 



and the density of states gitd) 



(8) 



The equation for self-consistency can be solved for large 
lattice sizes of the order of L = 30^ , that is indistinguish- 
able from the thermodynamic limit. 

In order to satisfy Ward identities, the calculation of 
conductivity within SCBA has to include ladder vertex 
corrections into the current-current correlation function 
of the Kubo formula.— However, for the model consid- 
ered in this study, these corrections are not effective, so 
that the conductivity is simply given by the product of 
spectral functions at different energies corresponding to 
the bubble diagram.— 

The quantities calculated within SCBA are reliable for 
not large el-ph couplings (up to A = 0.3) and for not too 
high temperatures (up to T = 0.4^^). 



B. Coherent potential approximation 

This approximation reduces to dynamical mean field 
theory in the limit of low particle density.^ Therefore, 
it is exact for infinite dimensions. Consequently, CPA 
self-energy is independent of the wave vector. 

CPA does rely on the entire displacement distribution 
and consists of the following self-consistent equation for 
the self-energy 



(9) 



A. Self-consistent Born approximation 

This approximation is valid in the weak el-ph coupling 
regime.^ It does not rely on the entire distribution of 
, but only on the first two moments: < >= and 
^ ^R^R i? "^ith r = ATtzX. Actually, tem- 

perature and el-ph coupling are cooperative to increase 
the interaction between electron and phonons. 

SCBA is an approximation for the averaged retarded 
Green function or equivalently for the retarded self- 
energy. For the considered model, the self-energy is in- 
dependent of the wave vector k and it is given by 



(6) 



where the Dyson equation provides the link between 
Green function and self-energy. Once the equation in fre- 
quency is solved, one can easily get the spectral function 
A{k, oj) through the following equation 



where the local Green function G^S*^ -j (lu) is given by 

Ri jR-i 



(10) 



A{k,uj) = -2SG'^*=*(/c,cj), 



(7) 



with /i chemical potential and electronic band of the 
model 

= —2tz cos {kz) — 2tx cos (kx) — 2ty cos (fcy). (11) 

For the solution of the equations, sizes of the order of 
L = 26^ are employed, getting results very close to the 
thermodynamic limit. 

Using CPA spectral functions, conductivity can be cal- 
culated. In the limit of infinite dimensions, there are no 
vertex corrections in the the current-current correlation 
function, so that only the bubble diagram contributes to 
the calculation of conductivity 1^ Therefore, while spec- 
tral properties calculated within CPA are reliable for any 
strength of el-ph coupling A and temperature T, the op- 
tical and transport properties are not very accurate, es- 
pecially in the regime of very low particle density. 
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C. Exact diagonalization 

Within ED, static and dynamic quantities of the sys- 
tem can be exactly derived. 

The electron part of the partition function can be 

calculated as (^Zei[{xj^ }]^, where Zei[{x^ }] is the 
quantum partition function of the electron subsystem 
given a deformation configuration |. This quan- 
tity is calculated through ED. The method provides an 
approximation-free partition function of the model in the 
semiclassical limit. The only limitation is due to the 
computational time being controlled by the L x L ma- 
trix diagonalization. This constrains our analysis up to 
L ~ 12^. In order to reduce the size effect, periodic 
boundary conditions are used. 

Within the same framework, it is also possible to cal- 
culate the spectral function, the density of states, and 
the conductivity averaging the electronic properties at a 
given ion displacement configuration over the entire set 
of configurations.^^ For instance, in the case of conduc- 
tivity, for each configuration, one calculates 



TT (l - I 



Re[(Tp,p{uj;{xj^J)] 
Y,Pril-Ps)\{r\Jp\s)\^S{Er-Es 

r^s 



(12) 



with p = x,y,z. In Eq. the Fermi distribution 

function pr 



Pr 



exp[(3{Er - /i)] + 1 



(13) 



corresponds to the exact eigenvalue Er, and {r\Jp\s) are 
the matrix elements between the exact eigenstates \r > 
and |s > of the current operator Jp along the direction 



Cp, defined as 



(14) 



We notice that, in contrast with spectral properties, the 
temperature enters the calculation not only through the 
displacement distribution, but also directly for each con- 
figuration through the Fermi distributions Pr- The nu- 
merical calculation of the conductivity is able to include 
the vertex corrections discarded by previous approaches. 
The mobility along the direction p is defined as 

Re[appiuj)] {M'^pA^-A^r,})] 
p = hm — = lim 

uj^o+ ec ui^o+ ec 

(15) 

In order to calculate observables for any finite lattice, 
we have to replace the delta function with a Lorentzian: 



-6{Er ~ Es+oj) ^ 



1] 



{Er ~E,+ a;)2 + ,^2 ■ 



(16) 



This replacement does not give any problem for quanti- 
ties at finite frequency, as soon as uj > r]. In the case of 
mobility, this procedure has to be correctly implemented 
since there is the limit w ^ 0. The correct expression 
is obtained performing the limits in the following order: 
w 0, 77 and L n- 00. This limit procedure is 
easily achieved within SCBA and CPA, due to the large 
values of L. However, for ED, it is less trivial. A detailed 
analysis has allowed us to obtain values for Wmm and 
rjmin that assure the determination of the correct value 
of the mobility for any temperature. For the ED mo- 
bility presented in this paper, we used uimin = lO^'^iz, 
77™„=3*10-2i^,andL = 123. 

In the following, we will measure energies in units of 
tz- Next spectral, optical and transport properties are 
discussed. 
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FIG. 1. (Color online) Spectral function within CPA at T = 
0.3 as function of frequency (in units of tz) for different values 
of kz — 2nn/{L)^^^ (in units of 1/a), with n — 1, .., 13 from 
left to right. In the inset, the width at half height as function 
of kz (in units of 1/a) for different values of temperature. The 
system size is L = 26^. 



III. SPECTRAL PROPERTIES 

Concerning spectral properties, focus is on the el-ph 
intermediate regime (A = 0.8), that is relevant for many 
organic semiconductors.^ 

In Fig. [TJ the spectral function A{k,uj) within CPA is 
reported at T = 0.3 for kz ~ 27rn/(_L)^/"^ and kx — ky — 
0. However, wc note that ED and CPA spectral functions 
are almost indistinguishable when L is small. 

With increasing the value of kz , the shape of the spec- 
tral function changes a lot. It is more peaked at small 
values of kz, while it broadens at the edge zone. For this 
reason, we show the widths at half height in the inset of 
Fig. [T] The width increases at the edge zone up to a value 
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very close to r/tz = 4AT. This quantity is related to the 
width at half height of the probability distribution for 
the random variable eg , which gives rise to the energy 
fluctuations of the electron system. Therefore, with in- 
creasing the temperature, the quasi-particles become less 
defined. Actually, aX T — 0.3, the behavior is more on 
the border of a quasi-particle description since the width 
is on average nearly of the order of t^. From the anal- 
ysis of widths at different temperatures, it is clear that 
there is a crossover from low to high T behavior and the 
increase at a fixed point of the Brillouin zone is linear as 
function of temperature, in agreement with experimental 
data^ 
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FIG. 2. (Color online) The quasi-particle energy within CPA 
as function of fcz (in units of 1/a) for different temperatures 
at A = 0.8. In the inset, the effective mass ratio m* /m as 
function of fcz (in units of 1/a) for different temperatures at 
A = 0.8. The system size is L = 26^. 



Another interesting point is about the peak positions. 
In Fig. [21 we report the quasi-particle energy dispersion 
derived by the peak positions of the spectral functions at 
different temperatures. Actually, the main peak of the 
spectral function inherits the bare band dispersion. With 
increasing temperature, only the peaks at small kz show 
a sensible reduction, but these corresponds to narrower 
peaks. Therefore, the widest peaks point towards bare 
band energies. This trend is also found in the effective 
masses. They are presented in the inset of Fig. [21 Only 
a factor of 2 is gained at center and border zone. More- 
over, effective mass values are quite stable with varying 
the temperature. Indeed, the system is not characterized 
by any quasi-particle with heavy effective mass. The be- 
havior of the quasi-particle band and mass is different 
from previous proposals based on polaron mechanismJ^ 
Indeed, if polaron formation takes place, a narrowing of 
the quasi-particle band and a marked increase of the mass 
is expected with increasing temperature. Actually, this 
behavior is exponential in the el-ph coupling and phonon 



distribution.-- In polaronic approaches, the entire bare 
band is reduced with a factor that is equal for all points 
of the Brillouin zone, while in our approach there is a 
strong difference in the renormalization at center and 
border zone. Moreover, in our approach, the increase 
of the width of the quasi-particle peak is playing a major 
role. 
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FIG. 3. (Color online) Density of states within CPA as func- 
tion of frequency (in units of tz) for different temperatures at 
A = 0.8. The horizontal lines depict the chemical potential 
at c = 0.01, c = 0.001, and c = 0.0001. The vertical dashed 
line at the energy —^At^ corresponds to the minimum of the 
bare band dispersion at the center zone. In the inset, the tail 
energy Etau as function of temperature at the same value of 
A. The system size is L = 26^. 

These behaviors are in agreement with ARPES data, 
that, as mentioned in the introduction, reflect the bare 
electron dispersions. Only in pentacene a very weak nar- 
rowing of the band has been found. From the analysis 
of the spectral functions shown in the previous figures it 
emerges that caution should be applied in the interpreta- 
tion of the data since width is another important parame- 
ter. It can increase with changing the system parameters 
making the band-like description not fully correct. 

The crossover from low to high temperature behavior 
affects the density of states extracted from the spectral 
functions. The marked width of the spectral functions 
gives rise to densities of states with a low energy tail. 
In Fig. [21 the density of states g is shown for different 
temperatures at A = 0.8. As shown in logarithmic scale, 
the density of states has an exponential behavior: g{uj) oc 
^-u/Etaii^ At very low temperatures,^^ this tail is well 
known and it corresponds to localized states. Moreover, 
it gives rough indications for the energy position of the 
mobility edge that should be close to the shoulder at the 
beginning of the tail. We notice that the energies in the 
tail increase very rapidly with temperature. 

In the inset of Fig. [31 the energy Etaii is reported as 
function of T . It shows a quadratic behavior as function 
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FIG. 4. (Color online) Density of states within CPA as func- 
tion of frequency (in units of t^) for different values of A at 
T = 0.3. The horizontal lines depict the chemical potential at 
c = 0.01, c = 0.001, and c = 0.0001. The vertical dashed line 
at the energy —SAt^ corresponds to the minimum of the bare 
band dispersion at the center zone. In the inset, the tail en- 
ergy Etaii as function of A at the same value of temperature. 
The system size is L = 26^. 

of T. At high temperature, it becomes of the order of half 
tz- The temperature trend and the order of magnitude 
of Etaii represent important pieces of information since, 
recently, there has been an attempt to estimate the tail 
in the density of stateS )^^'^" so that this study provides 
some clues for the interpretation of experimental data. 

It is interesting to analyze the role played by the chem- 
ical potential fi with varying the temperature. In Fig. |31 
the horizontal lines correspond to the values of the chem- 
ical potential at c = 0.01, c = 0.001, and c = 0.0001. For 
c = 0.01 and low temperature, /i is outside the tail char- 
acterized by localized states. Instead, for c — 0.0001 and 
low temperature, fi enters the energy tail and will pen- 
etrate into it with increasing temperature. In any case, 
with increasing temperature, /i goes into the energy re- 
gion of the tail. 

Similar trends are obtained when the density of states 
is analyzed as function of A. In Fig. HI the density of 
states as function of A is reported at T = 0.3. We notice 
that the tail is very small for A = 0.4, and it increases 
very rapidly as function of A. The tail energy Etaii is 
plotted as function of A in the inset of Fig. 01 We realize 
that Etaii is a linear function of A. This behavior in A 
again can be useful for the interpretation of experimen- 
tal data. Moreover, it clarifies that, in the intermediate 
regime, A and T play slightly different roles in contrast 
with SCBA, where they are completely equivalent in en- 
hancing the interaction of electrons with phonons. 

From the analysis of Figs. |3]and|ll it is clear that, in 
the regime of low carrier doping investigated in this pa- 
per, the chemical potential shows a strong tendency to go 
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FIG. 5. (Color online) Optical conductivity azz as a function 
of frequency (in units of tz) for different temperatures at A = 
0.3 and c — 0.01. In the upper inset, the width F of the 
conductivity as a function of temperature at A = 0.3 and 
c — 0.01. In the lower inset, the conductivity in the limit of 
zero frequency as a function of l/T at A = 0.3 and c = 0.01. 
The lattice size is L = 30^. 

towards the energy region of the tail with increasing the 
temperature and the el-ph coupling. This is a robust fea- 
ture in the intermediate regime. Therefore, this analysis 
allows to reconcile the results provided by the ARPES 
data (band-like description) with the trend that charge 
carriers appear more localized with increasing tempera- 
ture. This behavior is very important for the discussion 
of optical and transport properties in the next section. 

IV. OPTICAL AND TRANSPORT PROPERTIES 

As mentioned in the introduction, measurements of the 
optical and transport conductivity are important tools to 
investigate the properties of organic semiconductors. For 
this reason we have calculated these quantities within the 
studied model. First, we have analyzed the weak el-ph 
coupling regime within SCBA trying to understand the 
results of a band-like description. Then, we have consid- 
ered the intermediate el-ph coupling regime. We show 
the optical conductivity and mobility only along z direc- 
tion, since the corresponding quantities along the axis x 
and y can be roughly reproduced taking into account the 
anisotropy factors t^/tz and ty/tz, respectively. 

In Fig. [5l we show the conductivity for different tem- 
peratures at A = 0.3 and c = 0.01. For these values of 
parameters, the optical conductivity follows a Drude be- 
havior: a oc T/{uj^ + r^), with F the width of the Drude 
peak. In the upper inset, F is plotted as function of tem- 
perature: it follows a quadratic law. Finally, in the lower 
inset, we report the conductivity in the limit of omega 
to zero. In the adiabatic model that we are studying. 
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FIG. 6. (Color online) Optical conductivity a^z within ED as 
a function of frequency (in units of tz) for different tempera- 
tures at A = 0.8 and c = 0.01. In the inset, the conductivity 
a-zz in the limit of zero frequency within ED and CPA as a 
function of temperature at A = 0.8 and c = 0.01. Within ED, 
L = 10^ for conductivity at finite frequency, 12'^ for conduc- 
tivity in the limit of zero frequency. Within CPA, L = 26^. 



the temperature behavior is proportional to the inverse 
of temperature. This is the description of band transport 
in our model. 

The next step is to increase the el-ph coupling. The 
results from ED simulations are necessary to the anal- 
ysis since they clarify the effects of vertex corrections 
discarded by CPAi^ In Fig. [51 the optical conductiv- 
ity is reported for different temperatures at A = 0.8 and 
c = 0.01. First, we notice the reduction of conductivity 
in comparison with the data at A = 0.3. Moreover, the 
Drude peak does not longer represent the main spectral 
feature. Indeed, the conductivity exhibits a clear peak 
at low frequencies w < O.Giz, whose intensity decreases 
with the temperature moving slightly towards high ener- 
gies. This behavior can be ascribed to the role played by 
the vertex corrections into the calculation of conductiv- 
ity since the peak does not appear in CPA. Moreover, the 
peak positions do not depend on the broadening energy r\ 
(used in the calculation of the conductivity), that we still 
choose equal to 3 x 10~^t2. The result is of some interest 
since there are experimental evidence o^^'^*^ that, indeed, 
a peak is present at energies around Q2meV (500cm~^) 
lower than any charge transfer process.— This analysis 
confirms the findings of a previous paper for a related 
model.— 

In the inset of Fig. [6l we report the comparison be- 
tween the conductivity in the limit of zero frequency 
within ED and that within CPA for different tempera- 
tures. It is clear that, in the low temperature regime, 
CPA overestimates the value of conductivity pointing 
again towards the central role played by vertex correc- 
tions. The temperature behavior of Uzz is proportional to 



1/T, like in SCBA. At high temperatures, the CPA and 
ED conductivities at zero frequency tend to assume close 
values even if the finite frequency responses are quite dif- 
ferent. Furthermore, the crossover from low to high tem- 
perature behavior has a signature in the change of slope 
for both CPA and ED conductivities. The temperature 
behavior of conductivity appears to be proportional to 
1/T also at high temperature. 

Finally, we analyze the mobility fip along z direction as 
function of temperature with changing the particle den- 
sity. As shown in Fig. [71 the behavior is quite complex 
since it strongly depends also on the carrier density. For 
densities around one per cent, as pointed out in the previ- 
ous section, the chemical potential at low temperature is 
above the mobility edge. Thus, the mobility is metallic- 
like. However, the behavior at low and high temperature 
is different although all the range is characterized by a 
power law. Indeed, with increasing temperature, more 
and more states participate to the transport mechanism. 
However, the inclusion of localized states and itinerant 
damped states is only able to decrease the magnitude of 
the mobility. 

The situation is different for densities below one per 
cent. At those densities, the chemical potential is below 
the mobility edge, so that, at low temperature, we expect 
that the mobility shows an insulating character. With 
increasing the temperature, as emphasized in the previ- 
ous section, the chemical potential is always in the tail 
regions, but more states at higher energies (of itinerant 
nature) get involved providing the main contribution to 
mobility. Therefore, at intermediate temperatures, the 
mobility changes its character and exhibits a metallic 
power-law behavior. At higher temperatures, the mo- 
bilities look very similar for all the densities considered 
in this study. 

For temperatures T ~ wq, decreasing the densities, we 
observe that the chemical potential and, therefore, the 
most important states involved in the transport proper- 
ties cross the temperature dependent mobility edge (see 
Fig. 3). Translated in a polaronic framework, this would 
mean that the single particle sees a potential well due to 
el-ph coupling, hence, at very low densities, a polaron- 
like activated mechanism could come out. On the other 
hands, at higher temperatures, the scenario changes. The 
large broadening observed in the spectral function peaks 
at any density signals a new regime. As a result, the mo- 
bility decreases with the temperature acquiring a power- 
law behavior that do not depend significatively on den- 
sity. 

Finally, a comment on the behavior at very low tem- 
peratures (T << ujq) is in order. Indeed, in this regime, 
the adiabatic approach adopted in this work is no longer 
appropriate. As discussed in section II "The Model", 
at very low temperatures, quantum lattice fiuctuations 
could play a role. On the other hand, at the temper- 
atures T > Wo — 120K considered in this paper, our 
approach is fully consistent. This discussion clarifies in 
which sense the comparison with experimental data in 



8 



OFETs of rubrene and pentacene at very low densities 
and temperatures can be pursued. Of course, if the sam- 
ples are not very pure, but contain impurities or disloca- 
tions, then other mechanisms, not included in the present 
work, should be taken into account. However, our anal- 
ysis suggests that there is a narrow temperature window 
where an insulating behavior can be obtained without 
invoking the presence of impurities. In our view, the two 
different mechanisms could work together. 

Finally, the mobilities shown in Fig. [7] are expressed 
in our natural units, that is in terms of /xq = ea? /h. 
Taking a = ik^ one gets /uo ~ lcm'^/{V ■ s). At 
room temperature T = 0.25i2, the mobility is about 
3/io — 21cm? /{V • s), a value that recovers the right or- 
der of magnitude of experimental data in rubrene and 
pentacene.^ Moreover, the mobility results are consis- 
tent with those obtained within the picture of polaron 
formation. Actually, the power-law behavior as func- 
tion of temperature seems to be common to different 
mechanisms. In order to clarify the nature of the quasi- 
particles involved in the transport mechanism, it could 
be interesting to compare not only the mobility data, but 
also the optical and spectral properties of the present ap- 
proach with those that could be calculated in the frame- 
work of the polaronic approaches. Actually, at the mo- 
ment, these quantities in the polaronic framework are 
not available in the regime of parameters appropriate for 
organic semiconductors. Clearly, all the calculated prop- 
erties have to be discussed in relation with experiments. 
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FIG. 7. (Color online) Mobility within ED as a function of 
temperature for different particle densities at A = 0.8. The 
lattice size is L = 12^. 



V. CONCLUSIONS AND DISCUSSIONS 

Spectral, optical and transport properties have 
been calculated within the adiabatic three-dimensional 



anisotropic Holstein model recently proposed for the 
study of organic semiconductors. Different approaches 
have been used: SCBA, CPA, and ED. 

With increasing temperature, the width of the spec- 
tral functions A{k,uj) increases making the quasi-particle 
approximation less accurate. However, their peak posi- 
tions are not strongly renormalized in comparison with 
the bare ones. These results have been discussed in rela- 
tion with ARPES data that suggest a very simple quasi- 
particle energy dispersion. The marked width of the 
spectral functions is related to densities of states with a 
low energy exponential tail increasing with temperature. 
The order of magnitude and the behavior of the tail en- 
ergy as function of temperature and el-ph coupling have 
been analyzed in connection with recent experiments. 

The crossover between low and high temperature be- 
havior seen in the spectral properties affects also the op- 
tical and transport properties. In the weak el-ph cou- 
pling coherent regime, the optical conductivity follows a 
Drude behavior and the mobility can be interpreted in 
term of band transport. In the intermediate el-ph cou- 
pling regime, at low temperature, vertex corrections in- 
troduced through ED give rise to a temperature depen- 
dent low frequency peak in the conductivity, while, at 
higher temperatures, this peak broadens over large fre- 
quencies. This spectral feature agrees with the results of 
a previous paper for a related model;^ and compares well 
with recent experimental resultsJ^^— Finally, for high 
temperatures, the mobility is characterized by a power- 
law behavior as function of T. At low temperature, with 
decreasing the particle density, the mobility shows a tran- 
sition from metallic to insulating behavior. Results have 
been discussed in connection with available experimental 
data. 

In the model investigated in this paper, the mobility 
scales as the inverse of temperature. In experimental 
data of rubrene and pentacene, the mobility often de- 
creases faster as function of temperature. In order to 
simulate more realistic systems, the next step would be 
to consider Holstein-Peierls models?^"— where, as men- 
tioned in the Introduction, the lattice distortions affect 
also off-diagonal matrix elements of the electronic sys- 
tem. Actually, as shown also in recent papers by some 
of usj^i^ the non-local el-ph couplings in organic semi- 
conductors are not negligible, and they can affect the 
temperature dependence of carrier mobility. 

Long-range el-ph couplings could be included within 
appropriate models2^ in order to analyze the role of inter- 
face phonons in OFETsi^ Finally, it would be interesting 
to investigate the role of quantum lattice fluctuations. 
Their effects are small in the adiabatic limit, however, 
they could be important in the regime of very low tem- 
peratures. Work in this direction is in progress. 
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